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Recent cosmic ray, gamma ray, and microwave signals observed by Fermi, PAMELA, and WMAP 
indicate an unexpected primary source of e + e~ at 10-1000 GeV. We fit these data to "standard 
backgrounds" plus a new source, assumed to be a separable function of position and energy. For 
the spatial part, we consider three cases: annihilating dark matter, decaying dark matter, and 
pulsars. In each case, we use GALPROP to inject energy in log-spaced energy bins and compute 
the expected cosmic-ray and photon signals for each bin. We then fit a linear combination of energy 
bins, plus backgrounds, to the data. We use a non-parametric fit, with no prior constraints on the 
spectrum except smoothness and non-negativity. In addition, we consider arbitrary modifications to 
the energy spectrum of the "ordinary" primary source function, fixing its spatial part, finding this 
alone to be inadequate to explain the PAMELA or WMAP signals. We explore variations in the 
fits due to choice of magnetic field, primary electron injection index, spatial profiles, propagation 
parameters, and fit regularization method. Dark matter annihilation fits well, where our fit finds a 
mass of ~1 TeV and a boost factor times energy fraction of ~70. While it is possible for dark matter 
decay and pulsars to fit the data, unconventionally high magnetic fields and radiation densities are 
required near the Galactic Center to counter the relative shallowness of the assumed spatial profiles. 
We also fit to linear combinations of these three scenarios, though the fit is much less constrained. 



I. INTRODUCTION 

Several apparent anomalies in recent astrophysical 
data hint at a new source of high energy electrons, 
positrons, and possibly gamma rays, at the 10 GeV to 
1 TeV scale. The cosmic ray signals observed by Fermi 
[3-13 and PAMELA [| are direct evidence for these en- 
ergetic electrons and positrons (e + e~), which would lose 
their energy primarily through synchrotron radiation and 
inverse Compton scattering (IC). If the number density 
of these e + e~ rises towards the Galactic Center (GC), 
then this synchrotron and IC could explain the WMAP 
microwave "haze" @ and the Fermi diffuse gamma ray 
"haze" Q, respectively. 

It is difficult to explain these signals within the con- 
ventional diffusive propagation model and with standard 
assumptions about the interstellar medium (ISM). In this 
framework, the positron signal arises from secondary pro- 
duction from spallation of proton cosmic rays on the ISM. 
Assuming that 1. positrons and electrons have the same 
energy losses, 2. primary electrons and protons have the 
same production spectrum, and 3. the proton escape 
time decreases with energy, then the predicted positron 
fraction generically falls with energy, in contrast to the 
rising fraction observed by PAMELA. Katz et al. [8] 
point out these assumptions can be wrong, and explore 
alternative scenarios. Indeed, secondary production at 
shock fronts could explain the e + excess 0, E3] i but this 
would also imply an excess of anti-protons, which is not 
observed. We will not consider these alternatives further. 
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We examine here whether a new primary source of 
e + e~ is a viable explanation of the signals. First, the 
rise in the positron fraction measured by PAMELA sug- 
gests the presence of a new hard source of positrons [LL| . 
Second, the WMAP "haze" is consistent with a hard 
synchrotron signal in the inner galaxy, in addition to a 
soft-spectrum synchrotron component traced by Haslam. 
Though this decomposition is not unique, it is a good fit 
to the WMAP data. Third, the Fermi gamma ray "haze" 
similarly extends to \b\ > 30° above and below the plane 
in the inner galaxy. Neither haze correlates with the mor- 
phology of any known astrophysical objects or the ISM. 
(See Fig. Q]) 

Many attempts to explain the data operate by includ- 
ing a new component of high energy particles and gamma 
rays originating from one of the following sources: 

1. Annihilation of TeV-scale dark matter, 

2. Decay of TeV-scale dark matter, or 

3. An astrophysical source such as pulsars. 

These sources can produce energetic electrons, positrons, 
and gamma rays. In addition, the dark matter distribu- 
tion in the Galaxy is expected to be roughly spherical, 
providing at least qualitative agreement with the mor- 
phology of the gamma-ray and microwave haze. Never- 
theless, each explanation above has drawbacks. 

While annihilating dark matter may seem natural 
given a weak-scale WIMP which has a thermal freeze- 
out annihilation cross section, this vanilla scenario cannot 
explain the observed signals. Boost factors in the annihi- 
lation rate, arising from substructure or particle physics 
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0.097 


0.43 
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0.0765 


0.46 


15 


MED 




0.0112 


0.70 
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0.00595 


0.55 
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TABLE I. Typical propagation parameters consistent with 
low-energy cosmic ray data [32|. We use the "Default" pa- 
rameters and show the effect of using Ml and MED in Fig [5] 



enhancement, of order 100-1000 are typically needed, de- 
pending on the annihilation channels and dark matter 
mass. Significant model-building effort is also required 
to explain the lack of excess in the observed p/p flux 
[lH . For examples, see pl4l7j. 

In the decaying dark matter scenario, dark matter has 
the freeze-out annihilation cross section but also decays 
with lifetime t x ~ 10 26 s. These models also must explain 
why there is no excess in p/p, though no boost factors 
are required. Examples include [18l - [24l ] . 

The pulsar explanation is the least exotic, but there 
are significant astrophysical uncertainties in pulsar dis- 
tributions and e + e~ emission spectra. The Fermi cosmic 
ray signals can be explained by the presence of one or 
more nearby pulsars with hard e + e~ emission spectra 
[25U29I ] . However, pulsars are generally expected to be 
concentrated in the disk and it can be difficult to explain 
the shape of the WMAP and Fermi "haze" signals, which 
are much more spherical. See also [3(| EJ . 

In this paper we quantify how well each of these three 
scenarios can explain the data described above with- 
out resorting to model-dependent details of the particle 
physics or pulsars. Rather we use the data to determine 
the best-fit injection spectrum of electrons and positrons 
produced by each new source. We also show that a simple 
modification to the background electron injection can ex- 
plain the Fermi e + + e~ spectrum and the Fermi gamma 
ray spectrum but not the rest of the data. 

The standard procedure to analyse whether a model 
can explain the astrophysical signals is: 

pulsar or particle physics model 
spectrum of particles produced by the source 
propagation (e.g., GALPROP) 

comparison with data 

Often, one fits a specific dark matter or pulsar model 
to a subset of the "anomalous" signals described above. 
For dark matter, the particle physics model is usually 
processed through Pythia [33| to generate a spectrum of 
e + e~. The injection spectrum is the spectrum of e + e~ 
produced per unit source times the rate of production 
of e + e~ per source and the spatial distribution of the 
source. These injected e + e~ are propagated through the 
Galaxy to obtain a steady-state solution. The signals are 



then compared with data. 

While some analyses have studied the cases above in 
a less model-dependent way, the injection spectrum is 
assumed to have one of a few common forms (m , HU, [3|| . 

In this paper we effectively reverse the arrows in the 
procedure above. We fit the data from Fermi, PAMELA, 
and WMAP to expected backgrounds plus a new source 
which produces positrons and electrons. We assume the 
injection of the new source is separable in position and 
space. Rather than specifying the spectrum of e + e _ in- 
jected by the new source, we fit for this injection spec- 
trum as a function of energy. The shape of the spectrum 
is the same everywhere, and the spatial distribution is 
varied over several conventional models. Therefore, for 
our purposes, the three scenarios listed above differ only 
in their spatial distributions. 

In the pulsar case, the assumption that the pulsar in- 
jection is a separable function of position and energy is a 
crude approximation that allows us to fit the data with- 
out specifying the details of pulsar physics, since the po- 
sition dependence of pulsar populations and their e + e~ 
injection spectra is very uncertain. 

In our fits of the injection spectra, we simultaneously 
account for possible variations in the conventionally as- 
sumed spectrum and spatial distribution of the back- 
ground injected electrons, as well as propagation param- 
eters, magnetic field, and starlight densities. This takes 
into account the uncertainties in current models of the 
Galactic backgrounds. 

We describe the signals and their expected back- 
grounds in more detail in the next section. We then 
present the overall framework of the analysis. Predic- 
tions are computed using GALPROP, and we allow for 
variations in the background model. We then present 
the best-fit injection spectrum for each of the three sce- 
narios above, as well as the best-fit of the data to an 
arbitrary modification of the background electron injec- 
tion spectrum. Finally we present injection spectra for 
linear combinations of these scenarios. 



II. SIGNALS 

In this section we review the method of computing 
the signals and standard assumptions made in modeling 
the astrophysical backgrounds. However, in our fits we 
allow for variations in many of these assumptions. This 
is discussed in more detail in Sec. IIII Al 

In the conventional diffusive propagation model, the 
e~ cosmic ray density, dn(x,E)/dE, is the steady-state 
solution to the diffusion and energy loss equation: 

d Un{x,E Y 
dt V dE 

d_ 

dE 



= 



(1) 



, , 2 dn 



Q{E,x) 



where the first term represents diffusion, the second term 
energy loss, and the third term the source term. K(E) 
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FIG. 1. Maps of the (a) WMAP haze at 23 GHz and (b) Fermi gamma-ray haze at 5-10 GeV for the region -90 < I < 90 
and —45 < b < 45. The gamma-ray haze is obtained by subtracting the Fermi diffuse model from the data. All maps are 
centered on the GC. The data are compared to normalized maps of (c) pulsar synchrotron at 23 GHz and (d) synchrotron 
at 23 GHz from dark matter annihilation with an Einasto profile. The magnetic field has the form of Eq. [S] with tb = 6.5 
kpc. The morphology of the haze signals more closely resembles the signals from dark matter than from pulsars. We show the 
corresponding results for tb = 4.5 kpc in (e) and (f). Choosing a steeper magnetic field can change the morphology, but this 
is not preferred by the Haslam data (see Fig. [2]). The region —5 < b < 5 is masked out because we only fit the region b < —5. 



is the diffusion coefficient and b(E, x) is the energy loss 
rate. This equation holds separately for positrons. For 
both electrons and positrons, diffusive re-acceleration 
and galactic convection are negligible above a few GeV. 
Those effects are often relevant for other cosmic rays, 
which are governed by similar equations. We use GAL- 
PROP v50p.l to solve for steady state cosmic ray densi- 
ties. For a review, see (36[. 

For electrons, the source term includes primary elec- 
trons produced by supernovae and secondary electrons 
produced by collisions of proton cosmic rays on the ISM. 
We denote these sources by Qo(E,x). The source term 
can also include any new source of electrons, Qi(E,x). 
For positrons, the source term includes only secondary 
positrons and any new source of positrons. The spec- 
trum of injected secondary e + e~ is determined by the 
astrophysics of proton cosmic rays and their interactions. 

The injected primary electron spectrum is usually as- 



sumed to have the following energy dependence: 
(e 1m ,£<4GeV 

dN 

— oc ^ £> ,4 GeV < E < 2.2 TeV (2) 
[E 33 ,£>2.2TeV 

where j e can vary. dN/dE is the spectrum of e~ per unit 
source and is continuous. Eq. [2] has often been adopted 
in the past because the resulting cosmic ray fluxes ap- 
proximately agreed with the available data. Though we 
use this form as a default, we will also fit for an arbitrary 
modification to dN/dE. 

The number density for the supernovae that inject 
these electrons is commonly parametrized as 

n s (x) oc r a exp ( -/3— — ] @(r max - r) (3) 

V r Q .2 kpc/ 

where r is distance to the center of Galaxy, projected on 
the galactic plane, and z is distance perpendicular to the 
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galactic plane. The default GALPROP parameters are 
a = 2.35,^ = 5.56283, and r max = 15 kpc jaES]. 

The default normalization of the product n s x dNj dE 
is fixed such that the observed local flux from the primary 
electrons satisfies 

c din 

- — (34.5 GeV,z = O,r = r ) = 

3.15922 x l(T 7 (cm 2 • sr • s GeV)- 1 (4) 

which is consistent with the flux observed by Fermi. 

The diffusion of the injected e + e~ is governed by the 
diffusion coefficient, K(E), and L, the escape distance 
out of the galactic plane. K(E) represents the random 
walk of a charged particle in tangled magnetic fields, and 
is approximated as constant in space. It is generally as- 
sumed that K{E) = K (E/GeV) s , where K and 8 are 
propagation parameters. In Table p] we give some com- 
monly used values of Kq, 8, and L. [32j, |39(. Our default 
model assumes K = 0.097 kpc 2 /Myr, 8 = 0.43, and 
L = 4 kpc, though we will vary these parameters later. 
This choice matches cosmic ray data for protons, the B/C 
cosmic ray ratio, and was used in [40l |. 

As injected electrons and positrons propagate, they 
lose energy. The energy loss rate b(E, x) includes en- 
ergy loss mechanisms. The path length for an electron 
or positron losing 1/e of its energy is typically given by 
~ J KE/b, which is ~ 1 kpc around 1 GeV and be- 
comes shorter for higher energies, at least until the Klcin- 
Nishina limit [4lj . The dominant mechanisms for energy 
loss are IC scattering and synchrotron, where b(E) oc E 2 . 
Bremsstrahlung (brem), for which b(E) oc E, is sub- 
dominant above ~ 1 GeV and is far more localized to 
the disk. For a new high energy component of e + e~ , we 
neglect bremsstrahlung. 

The IC rate depends strongly on the model for the in- 
terstellar radiation field. We use the default GALPROP 
model [HI as a baseline. Meanwhile, the synchrotron de- 
pends on the model for the magnetic field. We assume a 
standard parametrization of the field, 

\B\ = B e^(-^y w (-^y (5) 

r is the distance to the center of the Galaxy, projected 
on the galactic plane. Typical parameters are Bq = 5/xG, 
tb ^5 — 10 kpc, and zb — 2 kpc. For our default prop- 
agation parameters, the GALPROP synchrotron predic- 
tion at 408 MHz best matches the Haslam 408 MHz map 
if tb ~ 8.5 kpc; see Fig. [2] 

This parameterization is consistent with observations 
of the large-scale (ordered) magnetic fields at 1-10 kpc 
J43j . The random component of the magnetic field is as- 
sumed to be proportional to the ordered fields, with a 



proportionality factor of approximately one [44|. Thus 
Eq. [5] is sufficient for our purposes, since our fits are not 
sensitive to the detailed structure of the magnetic fields. 
We increase or decrease the average strength of the mag- 
netic fields in the Galactic Center region by decreasing 
or increasing tb- 

The solution to Eq. Q] is the steady-state cosmic ray 
density, which then determines the photon signals. The 
gamma ray flux includes decay of 7r°s produced in proton 
cosmic ray collisions with the gas in the ISM, IC scatter- 
ing of e ± on interstellar photons, and bremsstrahlung of 
colliding with the ISM. The gamma-ray power in a 
given direction scales as: 

P^o oc J n gas {s) n p (s) ds , (6) 
-Pic oc J n*(s) n e ±(s) ds , (7) 

Pbrem OC J n gas (s) n e ±{s) ds. (8) 

where s is the coordinate along the line of sight. The ^23 
GHz microwave flux off the Galactic plane is primarily 
synchrotron radiation of electrons and positrons 

Psynch OC J \B{ S )\ 2 n e ±( S ) ds (9) 

where B is the magnetic field. 

A new source such as dark matter or pulsars can in- 
ject high energy electrons and positrons at 10 GeV to 1 
TeV. These new sources are included in Qi(E, x). In this 
paper we solve for the e + e~ injection spectrum which, af- 
ter propagation, yields the observed cosmic ray spectrum 
and gives rise to gamma rays and synchrotron radiation. 
Our fit will essentially determine Qi(E,xq), where x$ is 
the Earth's location. The spatial dependence of Q is fixed 
to be one of a few conventional models. 

These sources can also directly inject photons. There 
are primary photons from pulsars which are important 
at lower gamma-ray energies. Given our energy range of 
interest, we do not consider these further. 

In the case of dark matter annihilation or decay, gen- 
erally there are many channels through which dark mat- 
ter produces Standard Model particles, which can then 
decay on short time scales. The end products are e^, 
neutrinos, and photons. However, we do not consider 
these direct gamma rays any further. These can be pro- 
duced from 7r°s, final state radiation^ from r ± s or /i ± s, 
or a direct photon channel. For TeV-scale dark matter, 
these gamma rays can have higher energies than those 
observed by Fermi. Furthermore, in the fits below it is 
not difficult to produce enough gamma ray signal above 
10-100 GeV. In fact, direct gamma ray production can 
be rather constrained by observations [45|, EH • 



1 The documentation for GALPROP v50p incorrectly states that 2 In some papers, final state radiation is referred to as internal 
their parameter Bo is the magnetic field in the center of the bremsstrahlung. We use "bremsstrahlung" exclusively to mean 

galaxy. e + e _ cosmic rays colliding with the ISM. 
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Haslam 408MHz r B =4.5 kpc 




FIG. 2. (a) Haslam 408 MHz map. In the other panels we subtract the default GALPROP model for (b) tb = 4.5 kpc, (c) 
tb = 6.5 kpc, and (d) tb = 8.5 kpc from the Haslam map. The GALPROP model is normalized such that the total emission 
in the region £ £ [—10, 10], 6 € [—90, —5] matches the Haslam 408 MHz intensity in the same region. The constant offset is 
subtracted from the Haslam 408 MHz data. Note that local features like the North polar spur are not modeled in GALPROP 
and hence are not fit. 



A. Data 

We fit to the following data: 

• PAMELA J(e+)/(J(e~) + J(e+)) positron fraction, 
which displays a steep rise from 10-100 GeV [j| 

• Fermi (e + +e~) cosmic ray spectrum, which shows 
a slight hardening of the spectrum at a few hundred 
GeV [IHi 

• Fermi gamma ray spectrum, which shows a hard- 
ening of the spectrum at around 10-100 GeV, av- 
eraged over the haze region —15 < £ < 15 and 
10 < |6| < 30. Note the pion signal has been sub- 
tracted from the data Our background models 
match the pion component, shown in Fig. 11 of Q- 
This is not affected by the inclusion of new sources 
of electrons. 

• WMAP synchrotron at 23, 33, and 41 GHz aver- 
aged over -10 < t < 10 for -90 < b < -5, in 
2 degree bins. We also fit to the same data aver- 
aged over 10 < \£\ < 30, which we call the "high 
£" region of the WMAP data and is incorporated 
to include morphological information from the mi- 
crowave haze, [fj 

These data describe the "anomalous signals", which 
suggest the presence of a new source of electrons and 
positrons, and possibly gamma rays, at roughly 10-1000 
GeV. 



The Fermi LAT collaboration has provided a reference 
model for the diffuse emission 47], a detailed fit that 
includes a reference GALPROP model for IC and mod- 
els for a number of residual local features giving rise to 
bremsstrahlung at lower energies. Since we are not study- 
ing the detailed structure of the diffuse gamma rays and 
because IC and pions dominate at high energies, it is 
sufficient for us to use GALPROP to model the diffuse 
gamma ray emission in the haze region. 

We also do not attempt to fit the Fermi gamma-ray 
spectrum near the Galactic Center region nor the Fermi 
isotropic gamma rays. The signal near the GC suffers 
from large uncertainties in both the dark matter pro- 
file and the astrophysical backgrounds. The isotropic 
signal is extremely sensitive to the halo mass function. 
Some recent analyses have used these sets of data to con- 
strain dark matter explanations of cosmic-ray signals, for 
a variety of dark matter models and spatial distributions 
[H El Eil3|. In Fig. [3] we show that for the best fit 
spectra and spatial distributions in this paper, there is 
little or no tension between the predicted total signal and 
the data in the GC region. 

Finally, we also include the following data in order to 
help regularize the fit at lower energies and higher ener- 
gies: 

• AMS e + + er cosmic ray spectrum below 10 GeV 

m 

• HESS e + + e~ cosmic ray spectrum above 900 GeV 

m 
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FIG. 3. The fits in this paper are consistent with Fermi 
observations of the gamma ray spectrum in the region \£\ < 3, 
|6| < 3. Point sources have been subtracted. Solid colored 
lines show the predicted total signal for the best fits of the 
new sources considered in this paper, including backgrounds. 
The dashed lines show the contribution from only the new 
source. 

Fitting to the data from AMS ensures that the back- 
ground models are consistent with the low energy cosmic 
ray data. 

We include systematic errors in our analysis and treat 
them as statistical errors because we do not have the full 
covariance matrix. The energy calibration error of the 
Fermi data points is l^L, but rather than effectively 
increasing the error bars, we allow for freedom in the 
normalization of the background, discussed in the next 
section. The 15% energy calibration error has been in- 
cluded in the error bars used for HESS. 



B. Solar Modulation 

Though our focus is on high energy data, we ensure 
that our results are consistent with the low energy (be- 
low ~10 GeV) cosmic ray data from PAMELA and AMS. 
However, this data is extremely sensitive to the very lo- 
cal propagation in the heliosphere. To relate the data 
to GALPROP predictions for the local interstellar (LIS) 
spectrum outside the heliosphere, it is necessary to ap- 
ply solar modulation corrections to e + + e~ spectra. The 
solar modulation correction in the force-field approxima- 
tion is 

ME) Jlis{E + <S>) 
E 2 -m 2 e {E + <S>) 2 - m? e 1 ' 



where $ is the solar modulation parameter and J is the 
differential intensity dn/dE [53(. Because of the uncer- 
tainty in the force field approximation, we reduce the 
weight of the PAMELA and AMS data points below 10 
GeV, effectively multiplying error bars by a factor of 3. 
This is adequate to stabilize the fits at low energy. 

The solar modulation correction is applied to the GAL- 
PROP outputs. We also use the correction when convert- 
ing the positron fraction data of PAMELA into a positron 
flux data, using the AMS data on the intensity of e + +e~ . 
This will allow the fit to be linear below. 

However, these two data sets correspond to different 
parts of the solar cycle. We thus apply an inverse so- 
lar modulation correction to the AMS data to obtain 
the unmodulated positron intensity. Denoting the solar 
modulation correction by 5$, then the positron signal is 
obtained from 

JpAM(e + ) = ( ^f] ) x 
\J(e ) + J(e+)J PAM 

^ MM (C s [ J M S (e + + e-)]) (11) 

where $'s are solar modulation parameters. *&pam * s 
the solar modulation parameter for the PAMELA elec- 
trons, which we allow to be different from ^^> AM for 
the positrons. This approximately captures the charge 
dependence of the solar physics, visible in the time- 
dependent positron fraction at lower energies. In the 
above equation we applied S 1 *- to the total e + + eT 

P AM 

signal. Because the positrons are at most ~ 10% of the 
total flux, this approximation is justified. 



III. FITTING PROCEDURE 

We fit for the e + e~ injection spectrum that, when com- 
bined with a background model, best matches the cos- 
mic ray, gamma ray, and microwave observations. The 
steady-state e + e~ density is linear in the source function 
Q(E,x), so we take a Green's function approach in en- 
ergy space. The spatial dependence is fixed to be one of 
a few conventional models in each of the cases below. 

We inject delta functions of e + e~ at various energies 
and compute the signal from each delta function with 
GALPROP. Since GALPROP is discretized, in practice 
this amounts to propagating an appropriately normalized 
bin of energy. For each of these delta functions, GAL- 
PROP computes the steady-state e + e~ spectrum as well 
as maps of synchrotron and IC radiation at various ener- 
gies. We solve for the linear combination of these outputs 
that best matches the data. The best-fit injection spec- 
trum solution is simply the same linear combination of 
delta function injections (or in our case, energy bins). 

We inject e + e~ via the source term Q(E,x) in the 
propagation equation, Eq. Q] For dark matter annihi- 
lation, dark matter decay, and pulsars, the new source 
function Q\{E, x) of both positrons and electrons can be 
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written as 



Qi(E,x) = 



dN ftP ^ p x^ /e 

'" r -l £2L- £e 

X m x 2 
dN -1 

, dE 'p "P 



, ann 
, decay 
, pulsar. 



(12) 



Here dN/dE is the spectrum of electrons or positrons 
produced per unit source, normalized such that all the 
power per unit source goes into electrons!! p x (x) and 
m x are the energy density and mass of the dark matter. 
(<7t>)o is the thermal freeze-out cross section for annihi- 
lation, 3 • 10 _26 cm 3 /s. BF is a boost factor (from either 
particle physics or astrophysics such as substructure en- 
hancement). t x is the lifetime in the case of dark matter 
decay. t p and n p (x) are rate and density parameters as- 
sociated with pulsar emission rate and number density. 
Finally /g = /fi(e + + e~) is the fraction of energy going 
to electron-positron pairs. If /e = 1, then the total en- 
ergy of the electrons will be equal to m x for dark matter 
annihilation and m x /2 for dark matter decay. 

We also consider arbitrary modifications to the energy 
dependence of the background primary electron injection, 
Eq. [2] To accomplish this, we include an extra source of 
only electrons which has the same spatial distribution as 
the supernovae, Eq. [3J 



. dN 



, supernova 



(13) 



where r s is an arbitrary rate parameter that is fixed by 
matching to the data. 

Because we do not wish to a priori specify model pa- 
rameters, we instead implement the scenarios above with 
the following electron injection: 



i{E,x) = I 



Qi(E,x ) 
Qi{E,x ) 
Qi{E,x ) 
Qi(E,x ) 



Px( s o) , 

Px( S ) 
Px( s o) 

n p (x) 
n p (x ) 

n s (x) 
n s (x ) 



ann 



decay 
pulsar 
supernova. 



(14) 



where the local injection, Qi(E,xq) will be determined 
by the fit (xq is our location in the galaxy) . The positron 
injection is the same, except in the case of the source 
injection where there are no positrons injected. Only the 
spatial profiles distinguish dark matter annihilation, dark 
matter decay, or pulsars, in our fits. 

We bin the energies of the new source, Qi(E,x), and 
treat the particles in each energy bin independently. For 
example, we generally consider the energy range ~ 5- 
5000 GeV with 17 log spaced bins. The propagation of a 
given injection spectrum is just a linear combination of 
the propagation of each of the energy bins. 



3 The specific condition can be found in Sec. IIV Al for dark matter 
annihilation and Sec. IIV 51 for dark matter decay. 



The problem can be treated linearly because high en- 
ergy e + e~ are a tiny perturbation to the matter and ra- 
diation of the Galaxy. High energy e + e~ also almost 
never interact with each other or other cosmic rays; they 
dominantly interact with the ISM, radiation, and mag- 
netic fields. In GALPROP, the magnetic field is fixed 
and the usual feedback between cosmic rays and B field 
is absent. In this limit, the propagation of the individual 
energy bins is independent. 

We use this linearity to invert the propagation problem 
and determine the injected spectrum Qi(E,x ), given 
some assumptions about the spatial density of dark mat- 
ter or pulsars. Define the vector x by Xi = Qi(Ei,xo) 
for energy bin E^. The injection everywhere else is de- 
termined by the assumed spatial distribution. Also, let 
bj be the jth data point minus the galactic background, 
computed by GALPROP, for that data point. 

For each Xi, we propagate the injection and obtain 
a signal Aji. Thus A is a matrix which maps x to the 
predicted signal, and the columns of A give the predicted 
signal from each energy bin. We wish to compare the 
signal from the new source, A ■ x, with the background- 
subtracted data, b. 

The best fit x is determined by a goodness-of-fit test, 
which for a linear problem is a quadratic in the fit pa- 
rameters: 



X 



(A-x-b) T C- 1 (A-x-b) 



(15) 



where C is a covariance matrix containing the errors on 
the data. This is just a quadratic minimization problem. 
Note that we also include several other parameters in x 
that allows us to slightly modify the background predic- 
tions and improve the fit. This is described further in 
the following subsection. 

Finally, it is possible to obtain dN/dE and the other 
parameters in Eq. [12] This will be possible for the 
dark matter scenarios with additional constraints on 
/ E dN/dE dE and where dN/dE cuts off. This is dis- 
cussed more below when we describe the results for the 
scenarios above. 



A. Uncertainties 

The predictions for the signals discussed in this paper 
can have significant theoretical and astrophysical uncer- 
tainty. To capture the effects of these uncertainties, we 
include several parameters in the fits that essentially al- 
low for (small) variations in the background model. 

The main uncertainties are in the background primary 
and secondary e + e~ injection, since we are fitting all the 
data that could possibly constrain this. (For a more de- 
tailed discussion of uncertainties in the cosmic rays sig- 
nals, see dHHH.) 

In the primary injection, we allow 7 e in Eq. [2] to vary 
in discrete steps. We also fit for an arbitrary normaliza- 
tion factor N p relative to the condition in Eq. Usually 
we find N p ss 1.0 because the condition was chosen to 
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FIG. 4. All fit results for the three scenarios, over a 3 x 3 grid in background electron injection index (y e — 2.45, 2.55, 2.65) 
and scale for the magnetic field tb = 4.5, 6.5, and 8.5 kpc. These spectra were obtained from non-negative fits; the interpolated 
local injection density is plotted. Despite a wide range of assumptions about the background model, the results remain the 
same, qualitatively, for each scenario. 



approximately match the Fermi cosmic ray data. The 
source spatial distribution for primary electrons, Eq. |3l 
is also uncertain. Rather than considering the full range 
of possible source distributions, we simply allow for a dif- 
ferent normalization factor of the primary electron spec- 
trum near the Galactic Center, relative to Eq.[4] Because 
the diffusion length is ~ 1 kpc, this will not affect the lo- 
cal cosmic ray signal. This normalization factor, Nh, is 
fixed by matching the synchrotron background predic- 
tion onto the Haslam 408 MHz synchrotron map [56| for 
—90 < b < —5, averaged over —10 < I < 10. The contri- 
bution from the new high-energy source is negligible at 
this frequency. 

As for the secondary e + e~, it was shown in (32j 
that modifying the propagation parameters effectively 
changes the overall normalization of the local steady- 
state secondary positron flux by up to an order of mag- 
nitude. Thus rather than scanning over a large set of 
propagation parameters consistent with all the low en- 
ergy cosmic ray data, we allow the normalization of the 
secondaries to be a fit parameter, AT S Q 



There are several other adjustable parameters that can 
improve the fit and allow for variations in the background 
model. We give the complete list below. None of these 
will change the predictions for other cosmic ray data. 

The following parameters characterize the uncertain- 
ties of our theoretical models. We fit for these parame- 
ters simultaneously with the injection spectrum, as their 
effects can also be treated linearly: 

• Njc- The normalization of the background IC sig- 
nal near the center of the galaxy, relative to the 
GALPROP prediction. There are many uncertain- 
ties in the starlight density and spatial variations 
in the primary electron density near the galactic 
center. 

• AS": Zero-points of the WMAP signal. We allow a 
different one for each data set. 

• iV s : Normalization of secondary local electron spec- 
trum, relative to the GALPROP output for our 



4 We show the effect of changing propagation parameters on some 



fits in Fig. \E\ The result does not differ significantly relative to 
the error bars. 
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choice of propagation parameters. This can vary 
by up to an order of magnitude given theoretical 
uncertainties (55l. \El\. 

• N p : Normalization of local primary electron spec- 
trum, relative to Eq. [2 As mentioned above, this 
factor does not have to be the same as Nh- 

We include these parameters in x, and A is enlarged 
to include extra columns corresponding to each of these 
background signals. 

The signals are not linear in the following parameters, 
so we scan over a discrete set of these: 

• rs- The r-scale of the galactic B field, where the 
local B field is fixed to 5 [iG. See Eq. [SJ We include 
rg = 4.5,6.5, and 8.5 kpc, corresponding to B = 
33, 18, and 14 /iG in the center of the Galaxy. We 
used zb = 2 kpc. 

• 7 e : The index of the primary electron injec- 
tion spectrum above 4 GeV. We include j e — 
2.45, 2.50, 2.55, 2.60, 2.65, and 2.70. 

• ®ams,&pam'^pam : Solar modulation parame- 
ters for AMS and PAMELA, in the force-field ap- 
proximation, as described in Section [II B[ 

Though we allow these to be fit parameters, clearly in 
reality they have some definite form independent of our 
model. In Fig. 2] we show the best fit for a grid in tb 
and 7 e . Though the spectra do change, the qualitative 
features remain roughly the same. 

B. x minimization and regularization 

We are minimizing 

X^fA-x-bfC^fA-x-b) (16) 

where x is a vector of parameters we fit for, containing 
the injection spectrum as well as the normalization pa- 
rameters and offsets described above. C is a covariance 
matrix, so it is symmetric and positive. It can then be 
shown that the matrix A T C _1 A is positive semi-definite. 

Ideally the spectrum we derive is smooth and non- 
negative. However, the existence of null (or nearly null) 
eigenvalues of A T C _1 A means that there are directions 
in the parameter space where we can modify the spec- 
trum by large values with little change to the observed 
signals. This corresponds to, for example, changing the 
spectrum for two neighboring energy bins by a large pos- 
itive and negative amount respectively, such that the ob- 
served signal remains nearly the same. 

We regularize the spectrum by using only 17 log-spaced 
bins between 5 GeV and 5000 GeV. We also perform a 
non- negative quadratic fit following the algorithm in [58j j . 
All of the fit parameters should be positive except for A J 
which we find is always of the same sign for our data, so 
we can choose a convention where it is positive. 



To obtain errors on the spectrum, we find the eigen- 
vectors and eigenvalues of A T C~ 1 A. This allows us to 
change basis from the parameter space in x to a new pa- 
rameter space y where x 2 is separately parabolic in each 
parameter. The variance of these new parameters y is de- 
termined by computing the allowed shift of each param- 
eter, relative to the best fit, such that A^ 2 = ±1- Even 
though a non-negative constraint was imposed for the 
best fit, we consider the entire space of solutions within 
Ax 2 = ±1. 

Each linear combination of the so we sum 

the squares of the contribution from each iji to find the 
variance in Xi . The quoted error on each Xi is the square 
root of the variance. Because we are performing a non- 
negative fit, however, the positive and negative errors can 
be different. 

In Fig. [5] we show the entire range of possible varia- 
tions of the best fit injection spectrum with A^ 2 = ±1. 
We add to the best-fit spectrum all possible variations 
along the eigendirections, or all independent variations 
of yi. (We do not show the background normalization 
coefficients and WMAP offsets, though they are simulta- 
neously varying with the injection spectrum.) 

We also considered several alternative methods of reg- 
ularization, rather than non-negativity. As an example, 
we can impose smoothness by adding terms to \ 2 '- 

xlf f =X 2 + Vi (D E • x) T • (D E • x) 

+ % (D| • x) T • (D| • x) (17) 

where De and D E are finite difference and second- 
difference matrices, respectively. These matrices only 
act on the injection spectrum and not the other fit pa- 
rameters such as normalization factors and WMAP off- 
sets. 771 and 772 are tunable parameters that control the 
smoothness of the fit. In Fig. [5] we show the best fit us- 
ing this regulator instead of the non-negative regulator 
above. For an appropriate range of r\\ and 772 the so- 
lution is qualitatively similar to the non- negative result. 
Similarly, we tested several other regulation techniques, 
such as suppressing variations in nearly null eigendirec- 
tions. Again, for "reasonable" regulators, the result is 
qualitatively similar. 



IV. RESULTS 

We determined the best-fit injection spectrum for 350 
data points from Fermi, PAMELA, WMAP, AMS, and 
HESS. There are 29 fit parameters coming from 17 energy 
bins, 3 normalization factors, 6 WMAP offsets, and 3 
solar modulation parameters. Including tb and 7 e , then 
there are 31 fit parameters. Our results are summarized 
in TableGUand in Figs. 011 The details of the fit results 
for each scenario can be found in the following sections. 

In each of the following figures, we show the fits to the 

• e + e" flux data from Fermi, AMS, and HESS 
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FIG. 5. In the left column we show the space of possible solutions within A\ 2 = ±1 of the best-fit solution, which was obtained 
from a non-negative fit. The red lines are the best fits for (a) dark matter annihilation, (c) dark matter decay, and (e) pulsars. 
The spectra shown are interpolated between the bins (marked by red stars). In the right column we show other best fits 
obtained from using different propagation parameters, given in Tabled] or a different fitting regulator that enforces smoothness 
of the solution (from Eg. I17[) . Our results are robust to very different propagation parameters. The fit for M2 is not shown 
because L is only 1 kpc. Since the haze signals extend out to ~4 kpc or more, it is impossible for this set of propagation 
parameters to produce the haze. 



• positron flux obtained from combining the AMS 
and PAMELA data in Eq. [TP, 

• positron flux fraction J(e + )/(J(e~) + J(e + )) from 
PAMELA for comparison, though we did not di- 
rectly fit to this data 

• WMAP synchrotron data at 23 GHz and 33 GHz, 



and 23 GHz at high £; the data for 41 GHz and the 
high I data for 33 GHz and 41 GHz are included 
in the fit, but not shown because the fit looks ex- 
tremely similar to the plots already shown 

• Haslam 408 MHz data, used to fix Nh, as discussed 
in Sec. HHAl 
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TABLE II. Best fit parameters for annihilating dark matter, decaying dark matter, and pulsar cases to 350 data points. Annl 
and Ann2 had nearly the same x 2 but had different r_g so both results are displayed. In the supernova (SN) injection case 
there were 91 data points. We obtained mass, boost factor, and lifetime parameters from the best fit. In the last two columns 
we show fit results for linear combinations of these three scenarios. The fit errors on the normalization parameters N are less 
than 5-10% and thus are not shown. ^For the combination cases, we fixed the dark matter profiles to have Einasto a = 0.17. 



• Fermi gamma ray data, where the ir° background 
has been subtracted Q 

along with the best-fit local injection, E 2 Qi(E,xq). 

Before discussing the fits in detail, we emphasize that 
the results in Table [TT| and the spectra plotted here are not 
meant to be taken as precise answers but as qualitative 
guidelines for the necessary spectra, for each scenario, in 
order to explain the data. As shown in Fig. |4l the spectra 
vary with the background model, but the general features 
remain the same. Errors and variations in the solution 
were discussed in Section UTI Bl In addition, the effect of 
changing propagation parameters is shown in Fig. [5] 

Specific bumps and features in the spectra we find are 
more likely signs that the smooth background models we 
have assumed are not adequate. If there is any large 
systematic or unmodeled effect in the Fermi cosmic ray 
data, for example, it can change the features in best fit 
spectrum significantly. In particular, note that the shape 
of the high energy region of each spectrum above ^100 
GeV is only constrained by the high energy Fermi cosmic 
ray data since the Fermi gamma ray data is primarily 
only sensitive to the total power in this energy range. 
The other data are almost completely insensitive to such 
high energy particles. Thus the error bars on these bins 
are typically the largest. Furthermore, the high energy 
spectrum is more sensitive to changes in 7 e (see Fig. 0]). 

The low energy part of the spectrum is more severely 
constrained by all of the data. However, the low energy 
spectrum is also extremely sensitive to the bumps and 
features in the Fermi cosmic ray spectrum at low energies. 
This is very likely a sign that some features in the Fermi 



cosmic ray spectrum have not been properly included in 
the background model. For example, in Fig. [IT] we show 
a fit which allows both dark matter annihilation and an 
arbitrary modification to the energy dependence of the 
supernovae-injected electrons. The low-energy features 
can be fit by a modification of the supernova electron 
spectrum, while dark matter annihilation is still neces- 
sary to explain the signals above 10-20 GeV. 

Fig. [5] shows the "supernova" fit of the low energy 
PAMELA data, all of the e + e~ data, and the gamma 
ray data to a disk-like source with only electrons. This 
corresponds to a modification of the background primary 
electron spectrum and is implemented using the injection 
in Eq. 1141 The best fit spectrum we found is a hardening 
of the injection up to 1 TeV. Though this source modi- 
fication can easily match the cosmic ray or IC data, the 
disk-like spatial profile and lack of positrons produced 
are starkly inconsistent with the synchrotron signal and 
the PAMELA data. A new source is required. 



A. Annihilating Dark Matter Results 

The form of the injection for annihilating dark matter 
was given in Eq. [TT] and Eq. Q2J We assume the local 
dark matter density is po = 0.4 GeV/cm 3 (59j . 

Conventionally used dark matter halo density profiles 
are obtained by simulations and can be approximated 
by an Einasto profile, with 0.12 < a < 0.22 and a ~ 
0.17 on average [60|. This does not include substructure 
effects which can modify the effective spatial profile used 
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FIG. 6. Best fit to a modification of the primary electron injection spectrum, with Xred ~ -51- Black dashed lines are the 
background prediction for a model with j e = 2.65 and r_g = 8.5 kpc, though in this case we are fitting for the true background. 
Red dashed lines give the contribution of the new source injection, and solid red lines are the total. The gray shaded region is 
the error estimate on the best-fit injection spectrum. We have not attempted to fit the PAMELA data or the WMAP haze, 
which are difficult to produce. 



in Eq.H2l as in (6l|. 

We allow values of a = 0.12,0.17, and 0.22, with a 
core radius of r_2 = 25 kpc. In practice the shallower 
profile with a = 0.22 is always the best fit to avoid over- 
producing the gamma ray signal. These profiles only 
differ by a factor of ~ 2 at .1 kpc from the center of 
the galaxy. Though NFW profiles are also commonly 
used, their signatures can be approximated by one of 
these Einasto profile. We also considered spatial profiles 
which were Einasto squared times an r 1 / 4 or r 1 / 2 scal- 
ing, corresponding to an r dependent cross section [62j . 
Using these profiles can improve the \ 2 by 5-10, but the 
injection spectrum does not change significantly. 

In the annihilating case we found best fits with mag- 
netic fields of rs = 4.5, 6.5, and 8.5 kpc, all with 
X 2 around 140 and xled ~ -44- Conventional magnetic 



field models have rs closer to 8.5 kpc. Furthermore, in 
this case, the normalization factors N are ~1, so that 
the model is self-consistent. Thus we show the fit with 
rs = 8.5 kpc in Fig. [7] In Table [TT1 we give the fit param- 
eters for rs =6.5 kpc under the column "Ann2". 

For all three magnetic fields above, we found that an 
injection index of 7 e = 2.5 for the primary electron signal 
optimized the ratio between the PAMELA and the Fermi 
e + + er apparent dark matter components. However, for 
rs = 8.5 kpc the fit does not match the PAMELA data 
as well, as an excess of cosmic rays above 100 GeV can 
produce too many gamma rays through IC scattering. 

We can estimate several model parameters from the 
best-fit spectrum by relating Eq. [TJ] and Eq. [T2] To 
find the dark matter mass, we assume dN/dE cuts off at 
around m x . Though this estimate of m x depends on the 
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FIG. 7. Best fit for the annihilating dark matter scenario, with Xred ~ -44. The spatial profile of the dark matter is Einasto 
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lines give the contribution of the new source injection, and solid red lines are the total. The gray shaded region is the error 
estimate on the best-fit injection spectrum. 



rather uncertain high-energy part of the injection spec- 
trum, values of roughly 1 TeV are expected given the 
turnover in the e + + e~ data around 600-1000 GeV and 
the turnover in the gamma-ray spectrum at 100-200 GeV. 

Next, dN/dE was defined such that that the total en- 
ergy of the emitted particles sums to the mass of dark 
matter: 



E^dE 
dE 



(18) 



Therefore, integrating the local injection multiplied by 
energy gives 



E Q x (E,x )dE = (av} BF 



(A,) 2 f E (e + + e-) 



(19) 



Given an estimate of m x , we can therefore estimate BF x 
/s(e + + e~) in terms of the best fit local injection and 
known parameters. 



B. Decaying Dark Matter Results 

For the decaying dark matter case, we assume the same 
range of dark matter density profiles as in the annihilat- 
ing case. Again, in practice we will be limited to the case 
where a — 0.12. This time a steeper profile is required to 
produce sufficient synchrotron signal to fit the WMAP 
data. 

The model parameters can be determined from Eq. [T2l 
and Eq. [TJ] We assume dN/dE cuts off at around m x /2 
this time. Again, this cutoff is rather sensitive to the 
high-energy part of the spectrum, which has large error 
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FIG. 8. Best fit for the decaying dark matter scenario, with Xred ~ -41- The spatial profile of the dark matter is Einasto with 
a — 0.12. Black dashed lines are the background prediction for a model with 7 e = 2.6 and tb = 4.5 kpc. Red dashed lines give 
the contribution of the new source injection, and solid red lines are the total. The gray shaded region is the error estimate on 
the best-fit injection spectrum. 



bars, but values of > 2 TeV are expected given the data. 
By definition, dN/dE satisfies 

/ E ^ dE = m * A (20) 

Again, we integrate the local injection multiplied by en- 
ergy, giving 

[ w r> (t? - \ -i Po fE(e + + e~) 

I E Qx{E,x ) dE = t x — . (21) 

This allows us to determine the dark matter lifetime over 
the energy fraction. However, note that in many cases, 
dNj dE does not cut off in the energy ranges we consider 
and the spectrum is essentially unconstrained at higher 
energies. Then we only obtain bounds on the mass and 
lifetime. 



The best fit is shown in Fig. [8] There is no clear mass 
cutoff in the best-fit spectrum, so the mass of the particle 
can be from ~ 4 TeV to greater than 16 TeV. 

Because in the decaying scenario the injected power is 
proportional to p x and not p x , generally it is harder to 
generate enough synchrotron and IC signal. Both of these 
signals are in regions at least 5 degrees off of the galactic 
plane. The steeper dark matter profile with a = 0.12 is 
not enough to produce the signals. 

We found tb — 4.5 kpc can increase synchrotron near 
the center of the galaxy, but this gives a somewhat uncon- 
ventionally high value of the magnetic field in the GC, 
33^G. Fig. [5] shows that tb = 4.5 kpc also gives the 
poorest fit to the Haslam data, especially compared to 
tb = 8.5 kpc. In addition, a somewhat large injection of 
low energy electrons and positrons is required. However, 
for this large magnetic field, the IC signal drops. Thus 
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FIG. 9. Best fit for the pulsar scenario, with Xred ~ -46. The pulsar profile is given by Eq. I22I Black dashed lines are the 
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injection, and solid red lines are the total. The gray shaded region is the error estimate on the best-fit injection spectrum. 



the normalization Nic is rather large, Nic ~ 2.4. Even 
for fits with tb = 6.5 kpc, it was necessary for Nic ~ 2 
to obtain sufficient IC signal. This corresponds to rather 
high starlight density. 

While it is possible that the decaying dark matter 
can also produce gamma rays directly or through FSR, 
these signals are typically at higher energies, above 10- 
100 GeV. In this case, the large Nic factor for the back- 
ground IC signal indicates that there is a depletion of 
gamma rays at low energies, below 10 GeV. 

Though the decaying scenario nominally gives the best 
X 2 ~ 130 and xled ~ -41> the large normalization factors 
demand a more self-consistent modeling of backgrounds 
and uncertainties in order to fully justify the goodness of 
fit. 



C. Pulsar Results 



The range of types of pulsars, their spatial distribu- 
tions, and their e + e~ spectra is not well determined. 
As a crude model, we posit some spatial profile for the 
number density of pulsars and assume the spectrum of 
electrons and positrons has the same energy dependence 
everywhere. Generally pulsars are concentrated in the 
galactic disk, making it difficult to produce the haze. In 
Fig. [T] we compare the morphology of the synchrotron 
haze produced by pulsars to that produced by dark mat- 
ter. It is possible, however, that certain types of pulsars 
have a more spherical distribution (63| . 

We consider the following range of pulsar profiles, 
which span those typically used in the literature. (See 
[30l H3] for examples and references.) Assume a density 
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profile of the form 

( \A\ r ^ 
rip oc exp y—— j n p\P) ( 22 ) 

where p — \J x 2 + y 2 and the origin is at the center of 
the galaxy. We include profiles with z p = 0.08,0.5, and 
2 kpc. For the radial profile, 

^ KeXP (-l5l^)- (23) 

In practice, the three cases above look nearly identical 
because of diffusion. Another commonly used profile has 
n p ~ ( oexp(— p/4.0 kpc). However we do not consider 
this option further because the suppressed density near 
the center of the galaxy makes it even more difficult to 
produce the haze. 



Qualitatively, the pulsar results, Fig.|9l are rather sim- 
ilar to the decaying dark matter results, though the fits 
are even worse because of the disk-like rather than spher- 
ical profile. The best fit has tb = 4.5 kpc with signifi- 
cant low-energy injection and large normalization factor 
Nic of 2.6. Though it is possible that pulsars can pro- 
duce many low-energy gamma rays, it is unlikely these 
gamma rays can compensate for the background gamma- 
ray signal being 2 — 3 times too low. For rs = 6.5 kpc 
or 8.5 kpc, the pulsar scenario cannot produce sufficient 
synchrotron signal. 



D. Combination Results 

We fit for linear combinations of annihilation, decay, 
pulsars, and supernova injections, which not surprisingly 
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can provide better fits and alleviate the problems of each 
individual scenario. However, this extra freedom means 
that fits are much less constrained. For the fits presented 
here, errors in the spectra are much larger the spectra 
themselves and are not shown. 

To avoid the large normalization factors necessary for 
dark matter decay or pulsars to fit the data, we only 
consider linear combinations that include dark matter 
annihilation. Thus enough synchrotron can easily be pro- 
duced with lower magnetic fields. We also fix a to the 
middle value of a = 0.17. 

Generally, the power from dark matter annihilation at 
10-100 GeV produces sufficient synchrotron for a large 
fraction of the WMAP signal, while the power at 100- 
5000 GeV produces part of the gamma ray signal. De- 
caying dark matter or pulsars can then freely produce the 
rest of the local cosmic ray signal without significantly 



impacting the predictions from near the galactic center. 
Because decaying dark matter and pulsars do not give 
very different answers we don't show the fit to all three. 

Fig. [TU] is the fit to annihilation and pulsars. The fit 
to annihilation and decay is very similar, and thus not 
shown. Fig. [IT] is the fit to annihilation and a modifica- 
tion of the background primary electrons. This case is 
perhaps the most realistic of the possibilities, if a new 
source is allowed. There is an unusual spike at low en- 
ergies which just allows the primary electrons to more 
precisely match the features in the Fermi low energy cos- 
mic ray data, without overproducing any other signal. 
We have not shown a fit to a combination of pulsars and 
a modification of the primary electrons because, as men- 
tioned above, for both scenarios it was extremely difficult 
to produce the haze. A linear combination of the two 
would not alleviate this problem. 



18 



V. CONCLUSIONS 

We have thoroughly examined the annihilating dark 
matter, decaying dark matter, and pulsar explanations 
of the recent anomalous cosmic ray, gamma ray, and syn- 
chrotron signals. We investigated whether each scenario 
can fit all of these data simultaneously. Our analysis 
is independent of the particle physics or pulsar model 
details of each scenario and only depends on the spa- 
tial profiles and background models. We determined the 
necessary injection spectrum of electrons and positrons 
in each case in order to reproduce the data, including 
the effects of different background models, propagation 
models, and solar modulation. 

Though decaying dark matter is the best fit, the large 
normalization factors suggest that it will be difficult to 
find a fully self-consistent model with decaying dark mat- 
ter that can explain the data, without changing some as- 
pect of our model by a large amount. In particular, it 
may be necessary to find either a radiation field model 
with roughly twice as much starlight to produce enough 
IC, a rather large magnetic field of 33/iG in the Galac- 
tic center, enormous amounts of low energy electrons or 
gamma rays injected, a much steeper dark matter profile, 
or a combination of these. 

Pulsars give the worst fit; the disk-like profile makes 
it nearly impossible to produce both the gamma ray and 
synchrotron signals. Much like the decaying dark matter 
case, this suggests that dramatic re- assessments of back- 
grounds and models are necessary to find a self-consistent 
interpretation of the data. 

Annihilating dark matter, however, has self-consistent 
fits with conventional astrophysical background models. 
Though we had to choose a somewhat shallower dark 
matter halo profile with Einasto a = 0.22, it is still within 
the current range of profiles found by simulations. Fur- 
thermore, we can satisfy the gamma ray constraints from 



the GC. The boost factors are ~ 70//e, which at first 
seems much lower than the boost factors of ~ 1000 often 
used in the literature. Several factors enter in this dif- 
ference: our use of the updated po = .4 GeV/cm 3 rather 
than po = .3 GeV/cm 3 [59(, the relatively hard spec- 
trum allowed by the fit, and our assumption that the 
cutoff of the spectrum is m x . Given these factors, our 
result of ~ 70/ f E is typical of the models discussed in 
the introduction of this paper. However, the shape of 
the spectrum, combined with a lack of 7r° or p produc- 
tion, may still be difficult to achieve in current particle 
physics models of dark matter. 

Our results should be regarded as qualitative guide- 
lines to injection spectra. The specifics will necessarily 
change as both Galactic models and data are refined. 
The WMAP "haze" data will be superseded by data 
from Planck (65[, while data from Fermi and PAMELA 
will improve. In addition, cosmic ray data from AMS- 
02 (6(| may also soon be available. If the data does not 
change substantially, and if current models indeed de- 
scribe Galactic propagation and interactions, then the 
qualitative results of this paper will remain valid. 
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